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Abstract 

We measure interaction forces between pairs of charged PMMA colloidal particles suspended 
in a relatively low-polar medium (5 < e < 8) directly from the deviations of particle positions 
inside two time-shared optical traps. The particles are confined to optical point traps; one is 
held in a stationary trap and the other particle is brought closer in small steps while tracking 
the particle positions using confocal microscopy. From the observed particle positions inside the 
traps we calculate the interparticle forces using an ensemble- averaged particle displacement-force 
relationship. The force measurements are confirmed by independent measurements of the different 
parameters using electrophoresis and a scaling law for the liquid-solid phase transition. When 
increasing the salt concentration by exposing the sample to UV light, the force measurements 
agree well with the classical DLVO theory assuming a constant surface potential. On the other 
hand, when adding tetrabutylammonium chloride (TBAC) to vary the salt concentration, surface 
charge regulation seems to play an important role. 
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I. INTRODUCTION 



Understanding the interactions between (charged) particles is of great importance from 
a fundamental point of view to many natural, biological and industrial processes, which 
require control over the structure, stability and many other properties of a dispersion. [H [2] 
In addition, due to their mesoscopic size, colloids provide an ideal experimental system for 
the investigation of questions related to many particle statistical mechanics in- and out-of- 
equilibrium and related to both structural and dynamical properties of condensed matter. It 
is experimentally possible to visualize and follow their positions in real time using microscopy 
techniques. The particles can then be tracked to reconstruct their individual trajectories. 
From these, the interaction forces can be calculated. According to the classical Derjaguin- 
Landau-Verwey-Overbeek (DLVO) theory, the interaction between stable pairs of charged 
colloidal particles suspended in a dielectric medium is a sum of a repulsive electrostatic part 
and a generally shorter ranged attractive Van der Waals part.[3J H] However, when DLVO 
theory is quantitatively confronted with experiments, some refinements are required. [SHE] 
For example, the presence of adsorption of charged species at the colloidal surface and the 
complex charging mechanisms in nonaqueous media with low polarity, [HI [TO] make it unclear 
as to what extent the electrostatic interactions can be described with a volume fraction and 
phase independent sum of pair potentials. [TH [12] Measuring interparticle forces directly in 
these systems should help our understanding of the underlying physics. 

Different experimental techniques have been advanced to directly or indirectly measure 
forces acting between charged objects from a few femto-Newtons to sub-pico-Newtons. Total 
Internal Reflection Microscopy (TIRM) can measure the interaction forces between a mi- 
croscopic sphere and a flat surface. [13] The Magnetic Chaining Technique (MCT) has been 
used to directly probe the force-distance profile between magnetic colloidal particles. [HJ [15] 
Optical Tweezers (OT) provide a powerful tool to optically manipulate colloidal systems 
and have been widely used in recent years to investigate colloidal interactions. [TEH2"U] 
It is also possible to measure interaction forces by inverting the pair correlation func- 
tion g(r) in the case of weakly interacting systems using the Boltzmann distribution, 
U(r)/k B T = -\n [<7(r)].[ID[2Il[22] 

In the present paper, we use optical tweezers and a Nipkow scanning disk confocal mi- 
croscope to measure interaction forces between pairs of charged PMMA colloidal particles 
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suspended in a relatively low polar medium (cyclohexyl chloride (CHC), e = 7.6) at different 
ionic strengths, going from very low (purified solvent) to high ionic strength (added salt). 
Two particles are trapped, one is held in a stationary trap and the other particle is brought 
closer in small fixed steps. From the observed deviations of the particle positions inside the 
traps we calculate the interparticle forces. 

II. EXPERIMENTAL 

In our study, we used sterically stabilized and fluorescently labelled poly-methyl- 
methacrylate particles (PMMA) [23] with diameter of a = 1.4 /zm and size polydispersity of 
3%. The refractive index of CHC is n = 1.462 at A = 1064 nm, below that of the particles 
n p = 1.494, which allows optical tweezer experiments. In this low-polar solvent, charge dis- 
sociation still occurs spontaneously. [2"4] contrary to truly apolar media that require charge 
stabilizing surfactants. [25J Electrophoresis measurements in the dilute dispersion showed 
that the particles carried a positive charge of about +500e (with e the elementary charge). 
We studied systems at a volume fraction <p < 0.001, determined by particle tracking. Four 
different samples were explored. The first one consisted of a PMMA dispersion in purified 
CHC using a method described elsewhere. [26J In fact, the method offers a very convenient 
way to quickly reduce the conductivity (from ~ 1000 pS/cm to less than 20 pS/cm) and at 
the same time, alumina is often used as a desiccant to remove traces of water. Additionally, 
the purified solvent is stored with added molecular sieves (4 Angstrom, Acros Organics) that 
serve as an adsorbent. In the three other systems, the ion concentrations were increased 
by exposing the initial purified sample under an UV lamp at different exposure times. It 
is known that UV exposure facilitates the partial degradation of cyclohexyl halides into 
ions, [27] but it is difficult to quantify the salt concentration by conductivity measurements 
without knowing exactly which ions are generated. However, our measurements give an idea 
of how the different parameters in the same sample vary when exposed to UV light. The 
dried particles were initially mixed with the purified CHC and allowed to equilibrate for 48 
hours before measurements were performed. The sample was then gently remixed, placed 
in a glass capillary of 0.1 x 2.0 mm inner dimensions (VitroCom) and sealed with Norland 
68 UV glue. While curing the glue, the sample was covered with aluminium foil to prevent 
solvent degradation. 
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The measurement of particle interactions was performed using optical tweezers. [28j 129] 
We use a pair of optical tweezers formed by focusing a 1064 nm laser, using the same lens as 
used for imaging. We applied time-sharing using accousto-optic deflectors to generate two 
traps and to vary their positions. [30j Particle imaging was performed with a Nipkow-disk 
scanning confocal microscope (CSU10, Yokogawa) and recorded on a digital video camera 
(Evolution™ QEi) as in Ref.jHj. A 100 x 1.4 NA oil immersion objective (Leica PLAN APO) 
was used. The particles were dyed with Rhodamine and were excited with a Millennia V 
diode-pumped laser beam (A = 532 nm). The trapped particles were located in the plane 
at about 14 /im above the bottom of the capillary glass wall to avoid possible effects of the 
sample boundaries. One of the traps was brought closer to the other one in small steps of 
about 500 nm every 2-3 seconds. For every step, 1000 images of 80 x 21 pixels were recorded 
to sample the Brownian motion of the particles inside the traps. All images were processed 
to extract particle positions using home-made software based on methods similar to that 
described by Crocker and Grier. [2"§] 

Consider two trapped particles. One (at position r) is held in a stationary trap (at the 
origin) and the other particle (situated at distance R from the first particle) is brought 
closer in small steps (see Figjl]). The particle in the stationary trap will feel both the force 
F W eii(r) exerted by its optical trap (due to a well potential C/ we u(r)) and a force F(R) due to 
the presence of the other particle. The additional force F causes the equilibrium position of 
the particle in the stationary trap to shift by a small amount. As the mobile trap is moved 
closer to the stationary one, the strength of this force increases and the particle deviates 
further. The particle in the mobile trap is similarly pushed in the opposite direction. For 
sufficiently large forces, Brownian motion can cause one of the particles to escape its trap. 
Figure [T] shows the positions of two particles as the mobile trap approaches, until one of the 
particles escapes. At large distances, the pair interaction between the particles is negligible, 
but as the traps approach, each particle deviates from its well center due to the electrostatic 
repulsion. If the trap potential is known, the equilibrium position of the particle in the trap 
can be calculated from the Boltzmann distribution. If we assume that the force F exerted by 
the second particle does not vary significantly over the volume where the particle fluctuates; 
F(R) F(R ± dr), where dr is the typical fluctuation displacement of the particle inside 
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FIG. 1: (Top left) A sketch of two particles confined to two optical traps. (Top right) Confocal 
microscopy images of the two particles as the mobile trap (MT) approaches the stationary one 
(ST). The red contours on the particles are centered on the positions taken from the tracking code. 
(Bottom) Trajectories of the two trapped particles (1 pixel = 0.181/nn) projected onto the x-axis 
(chosen along the line connecting the two traps) for different frames. Inset: a zoom of the particle 
position deviations inside the stationary trap. 

the trap, the expected position is given by: 

/ r exp (-/3 {U wc n(r) - F ■ rjj dr 
{T) /exp(-/?{tW(r)-F.r})dr ' 1 J 

where the integrations are carried out over the volume of the trap. After (numerically) 
calculating the integral, this relation can be inverted to calculate the force on a particle 
from its measured position. For perfectly harmonic traps, F oc (r). 

To calculate the well potential we used the Mie-Debye representation given in Ref . [32} |33] . 
In these calculations, the parameters describing the laser trap configuration are the beam 
opening angle 9 and 7 the ratio of the objective focal length to the beam waist u . We used 
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6 = 64.245°, 7 — 1-21 to model our laser spot with the lOOx 1.4 NA lens with appropriate 
overfilling. [Ml I3"5] The spherical aberrations, due to the refractive index mismatch between 
the coverslip and the CHC, depend on the distance from the geometric focus to the glass 
surface (14 fim). To calculate the potential energy of the particle in the trap, the trapping 
force was integrated in the lateral direction at the trapping depth. The calculations show 
that for the parameters used in our experiments a distinct deviation from harmonic behavior, 
a relative increase in the stiffness at larger deviations, is expected. [361 For the measured 
forces with the largest deviations, this effect can be up to 20%. By using the trapping 
potential from calculations we rest reassured that we do not underestimate these relatively 
larger forces. The trapping forces and potential assumed to scale with the power of the trap, 
such that the well potential is U we n(r) = aU C3 x c (r). The prefactor a was measured from the 
standard deviation of the particle positions at large distances. 



III. RESULTS AND DISCUSSION 



From our measurements, we plot the particle deviation in the stationary well (r) ver- 
sus the particle separation (R) (for an example, see Figj2|. Initially, the particle exhibits 
Brownian motion inside the trap of a magnitude less than 0.1 fim, i.e. < 12% of the par- 
ticle diameter as demonstrated by the probability distribution of displacements plotted in 
the inset of the figure. From the standard deviation of this distribution we determine the 
proportionality factor a, which sets the well depth. We then calculate average deviations 
at different interparticle distances by binning the data shown in Fig|2j and convert these to 
forces using Eq. Q. Performing measurements with different laser powers (200 to 800 mW) 
gave the same results within our experimental resolution, which clearly demonstrates that 
possible light-induced particle interactions can be neglected. [T§J [3SJ EH] 

Our data of the pair-interaction forces are compared with the prediction of the linearized 
Derjaguin approximation [2] with constant surface potential of the form 

V ; 4A B 1 + exp(-Ka(R/a - 1)) ' K J 

where = e^o/k-BT is the dimensionless surface potential, Xb the Bjerrum length in CHC 
(~ 7.3 nm), k = \/47t\bC s the inverse Debye length assuming only monovalent ions of 
concentration C B , a the particle diameter and R the distance between the trapped particles. 
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FIG. 2: Measured particle deviations inside the stationary trap plotted against particle separa- 
tions for the sample with the highest salt concentration studied in this work. The particle starts 
feeling the approch of the mobile trap at a particle separation of about 2a. Inset: the probability 
distribution of displacements inside the stationary trap before any deviation and the fitted line is 
a Gaussian approximation to the experimental data. 

Figure [3] shows the results for the different sets of measurements as the concentration 
of salt is increased by exposing the initial purified sample under an UV lamp at different 
exposure times: 0, 1, 12, and 16 hours respectively. When increasing the salt concentration, 
the range of the interparticle forces decreased. The data are well fitted with the DLVO 
force given by Eq. ^ with a fixed surface potential <p — 5-3 (or ^ ~ 135 mV) and a 
Debye screening length decreasing from ~ 1.3 fim in a purified sample to k^ 1 < 0.2 
/im in the sample with the longest UV exposure time (highest salt concentration). The 
interaction force F(R), initially soft and extremely long-ranged, becomes shorter-ranged for 
higher salt concentration, approaching that of hard-sphere systems. When performing the 
same measurements in the same sample with a different speed of approach, the resulting 
interactions do not differ significantly. This suggests that hydrodynamic effects due to the 
movement of the traps are negligible. 

Since we do not follow the deviations in the z-direction but measure the projection of the 
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FIG. 3: Interaction forces between charged colloidal PMMA particles suspended in cyclohexyl 
chloride with different salt concentrations resulting from different UV exposure times: (a) initial 
purified sample, (b) 1 hour UV exposure, (c) 12 hours and (d) 16 hours. The lines are fits based 
on the DLVO theory (Eq. Q) with fixed surface potential = 135 mV, giving values for Debye 
screening lengths k~ 1 of 1.3, 0.87, 0.35 and 0.22 //m respectively. 

particle separation onto the xy plane, we slightly underestimate the distance between the 
two particles. In addition, the assumption that the force is approximatly constant over the 
integration volume in Eq. Q will lead to a small systematic errors in the measured forces. 
We performed Monte Carlo simulations to estimate these effects, simulating two particles in 
two optical traps with the interactions and well shape based on our experimental parameters. 
After analyzing the resulting data in the same way as the experimental data, the deviations 
in the obtained values for kg and \&o were at most 3% and 1%, respectively. 

To further test the reliability of this method, we conducted two more measurements using 
two other dispersions containing a controlled amount of salt tetrabutylammonium chloride 
(TBAC, Sigma- Aldrich). We independently estimated the Debye screening lengths from the 
scaling law for the liquid-solid phase transition in Yukawa systems, |40j and measured the 
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FIG. 4: Interaction forces between charged colloidal PMMA particles in samples I and II with 0.026 
and 0.26 fiM added TBAC salt respectively. The lines are fits based on the DLVO theory with 
constant surface potential (full lines) with <f> = 5.64 and 3.37, giving values for screening lengths 
kct of 1.86 and 2.58 respectively, and with constant surface charge (dashed lines) with Q = 4.56 
and 2.55, giving values for screening lengths kg of 1.6 and 2.06 respectively. 

particle surface charges by means of electrophoresis. We prepared a solution of TBAC in 
purified CHC, which we allowed to equilibrate for a week. We then filtered the saturated 
solvent (~260 /xM TBAC) and diluted it by adding an amount of purified CHC to prepare 
two solvents I and II containing 0.026 and 0.26 /iM of TBAC, respectively. The measured 
force-distance profiles of PMMA particles in these two salt added solvents are again well 
fitted with the DLVO force but with different surface potentials. The data are also well 
fitted with the constant surface charge formula [21 Hi] (see Fig|4]and Table [I] for the different 
parameters). When the two samples are left vertically, sedimentation of particles induced 
crystallization. The samples were imaged with a tilted confocal microscope that allows to 
scan at all heights. The average particle-particle separation d for the two systems near 
freezing were 5.3 and 4.1 /iin respectively (see Fig|5]). Using the scaling law proposed in 
Ref. [3D], the coupling parameter needed for crystallization: 
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FIG. 5: Confocal images of the two sedimenting samples I and II with 0.026 and 0.26 fiM added 
TBAC salt, respectively, below (a), near (b) and above (c) freezing. The average particle-particle 
separations are 5.3 and 4.1 jjm in I and II respectively. Using the coupling parameter needed for 
crystallization (Eqj3j), this gives an estimate for the screening lengths ko of about 2.3 and 2.8 
respectively. 

where U(d) is the potential energy at the typical particle-particle separation. Using surface 
potentials from our force fits, this gives values for kg of about 2.3 and 2.8 respectively, in the 
same order of magnitude as from our force measurements. In addition, we obtained indepen- 
dent measurements of the ionic strength from the conductivity of the particle-free solvents. 
A commercial conductivity meter (Scientifica 627) was used to measure the conductivity of 
the two solvents with the added salt. This yielded conductivities of 1600 and 9800 pS/cm 
respectively. Using Walden's rule,jl2] the corresponding screening lenghts are kg — 2.1 and 
5.5 respectively. The screening parameters obtained from conductivity measurements are on 
the same order as those obtained from our force measurements and EqJHJ although the value 
for sample II is slightly higher. By fixing kg = 5.5 for sample II, we were neither able to fit 
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Constant Surface Potential (Eq. (2 
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Constant Surface Charge (Ref.0 E]) 
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1.6 
2.06 


4.56 
2.55 


547 
237 


4.38 
2.43 


980 
510 



TABLE I: Measured parameters for the two samples I and II with 0.026 and 0.26 /jM added 
TBAC respectively: the screening length kg, the dimensionless surface potential <j) = efo/^B? 1 , 
the dimensionless surface charge Q = eZ* / Aitsk-QT n~ l , and the effective charge Z* . The effective 
charges quoted for the constant potential fits are for infinite particle separation. The last two 
columns are based on the electrophoresis measurements. 

our force measurement data assuming a constant surface potential nor a constant surface 
charge. Therefore, we believe that the estimated value kg = 2.8 from Eqj3]is more accurate 
in this case. 

Electrophoretic measurements on the same dilute samples were conducted by driving the 
particles in a dc-electric field (E « 1 — 3 V/mm) and measuring the electrophoretic mobility 
from particle tracking of the confocal images. Using the values of the Debye screening length 
k~ 1 from our force fits and the measured mobilities, the surface potentials were obtained 
using recent calculations for electrophoresis. [4~3] This gives values of e i = 4.5 and 2.42 for 
samples I and II respectively when using the screening lengths from the constant surface 
potential fits and, 4.38 and 2.43 for the case of constant surface charge. Finally, to translate 
the surface potential into a particle charge Z*, we used the empirical relationship proposed 
by Loeb et a/.[H] The results of the electrophoresis measurements are slightly lower than 
the measured force parameters with a constant surface potential but slightly higher than 
the parameters obtained under the assumption of a constant surface charge. These results 
suggest the importance of charge regulation in the case of PMMA dispersed in CHC with 
added TBAC.[II] 
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IV. CONCLUSION 



In summary, we have measured the interaction forces acting between charged PMMA col- 
loidal particles suspended in a relatively low dielectric medium directly from the deviations 
of particle positions inside an optical trap. The measured forces include all possible effects 
such as the influence of the surrounding particles. Confocal microscopy can be used to track 
the particles at relatively high volume fractions ~ 15 % where optical tweezers can still be 
operated. As a result, this method allows to measure particle interactions at higher volume 
fractions. Additionally, our measurements quantified for the first time the effect of exposing 
a PMMA/CHC dispersion to UV light. When increasing the salt concentration by exposing 
the sample to UV light, the force measurements agree well with the classical DLVO theory 
assuming a constant surface potential. On the other hand, when adding tetrabutylammo- 
nium chloride (TBAC) to vary the salt concentration, surface charge regulation seemed to 
play an important role. Future work will explore the effect of particle concentration and the 
nature of the interactions in dense systems. 
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